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Abstract. The Bogdanov-Takens system has at most one limit cycle and, in the 
parameter space, it exists between a Hopf and a saddle-loop bifurcation curves. The 
aim of this paper is to prove the Perko's conjectures about some analytic properties 
of the saddle-loop bifurcation curve. Moreover, we provide sharp piecewise algebraic 
upper and lower bounds for this curve. 

&. 

1. Introduction 

The Bogdanov-Takens system 



C/3 

Q 



x =y, 
y' = —n + by + x 2 + xy, 



has been introduced in [H [T71 Q2]. It provides a universal unfolding of a cusp point of 

codimension 2 and it is considered in many basic text books on bifurcation theory; see 

for instance [31 El [12] . Some global quantitative properties of its bifurcation diagram 

are not known. In 1992, Perko stated two conjectures about analytic properties of the 

saddle-loop bifurcation curve in the parameter space; see [13]. The aim of this work is 

>• ', to prove both conjectures. 

The interesting bifurcations only appear in the region n > 0, because in this case the 

system has two critical points, (±-y/n, 0), a saddle and a focus. Therefore it is natural 

to introduce a new positive parameter m = y/n. So, we will consider the following 

CO ! expression of the above system 

O 



J x =y 

y' = —m 2 + by + x 2 + xy, with m > 0. 



Before presenting our results, we recall the known properties about the bifurcation 
diagram of system ([!]). All the qualitative information of this diagram and part of the 
quantitative one are known; see [TOl [131 US US]. In particular, it is proved in [10] that 
this system has at most one limit cycle and that when it exists it is hyperbolic and 
unstable. This information, together with the fact that system ([I]) is a rotated family 
of vector fields with respect to b, allow to show that the limit cycle exists if and only 
if b*(m) < b < m, for an unknown function b*(m). This holds because fixing m and 
decreasing b, a unique unstable limit cycle borns via a Hopf bifurcation for b = m, 
increases diminishing b, and disappears in a saddle- loop connection for b = b*(m). The 
corresponding phase portraits are drawn in Figure [TJ and a sketch of its bifurcation 
diagram is given in Figure [21 
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(i) b < b*(m) 






(ii) b = b*(m) (iii) b*(m) < b < m 
FIGURE 1. Phase portraits of system (pTJ) 



fiv) b > m 



Some quantitative information about b*(m) is given by Perko in [13]: 

(i) It is an analytic function. 

(ii) It holds that max(-m,m — 1) < b*(m) < m. 

(iii) At m = 0, b*(m) = 5m/7 + 0(m 2 ). This term is computed by using the Melnikov 
method; see also [8]. 

As usual, we write f(m) = 0(m p ) or g(m) = o(m p ) at m = m GlU {oo} if 

/(m) 



lim 



if e 



or 



lim 



0. 



The lower bound given in item (ii) is improved in j9] applying the Bendixson-Dulac 
Theorem and proving that max(m/2,m — 1) < b*(m) < m. 




FIGURE 2. Sketch of the bifurcation diagram of system (TTJ. The open colored 
region is the one containing the limit cycle 

Item (iii) has been improved recently in [7] using a different approach, based on the 
construction of algebraic curves with a loop that is without contact for the flow of the 
system. The authors obtain that, at m = 0, 

,*/ x 5 72 2 30024 3 2352961656 4 5 . 

h*( m ) = -m J [ vri m 6 m 4 + 0(m b ). (2) 

v ; 7 2401 45294865 11108339166925 v ; v ' 

The Bogdanov-Takens system for parameters in a neighborhood of infinity is also 
studied in [2]. The aim of that work was to understand the presence of the limit cycle 
for the system in terms of slow-fast dynamics. No quantitative information about the 
shape of the curve b = b*(m) is given there. 

In [T3] a different, but equivalent, expression of system ([T]) is considered. Next con- 
jectures correspond to Perko's ones translated to ([T]). 
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Perko's Conjectures ([13]): Let b = b*(m) the function that corresponds to the 
saddle- loop bifurcation curve for system (pQ). Then, 

(X) for m large enough, b*(m) = m — 1 + 0(-^=), 

(XX) it holds that max(5m/7, m — 1) < b*(m). 

To facilitate the reading of this work, original Perko's formulation of above conjectures 
is recalled in Section [6j 

Both conjectures are immediate consequences of Theorems [1] and [2j Moreover, The- 
orem [2] significatively improves the global lower and upper bounds given above. 



Theorem 1. Form large enough, b*(m) 

Theorem 2. It holds that 
5m 



m 



V m I 



max 



m 



1 < b*(m) < min 



(5 + 



m)m 



37 
12' 



7 + ^m 



m 



25 
7m 



To prove our results we develop the method introduced in [7J, adapting it according 
to small or large values of m. The basic idea is as follows: for given positive values of 
b and m such that b < m, we want to know if b > b*(m) or b < b*(m); or equivalently 
to prove the existence or non-existence of the limit cycle. Due to the uniqueness and 
hyperbolicity of the limit cycle these two situations can be distinguished constructing 
negative or positively invariant regions, as it is shown in Figure [31 and employing the 
Poincare-Bendixson Theorem. 




b > b*(m 




b < b*(m 



FIGURE 3. Negative and positive invariant regions around the focus 



The most difficult part of this approach consists in constructing these negative or 
positive invariant regions delimited by loops. In [7] , the closed loop around the attract- 
ing point to prove ^2} is proposed to be the loop of an algebraic self-intersecting curve 
whose vertex is on the saddle point and whose branches approximate its separatrices. 
Here we use and develop this idea when m is small. One of the key points for proving 
that b = 5m/ 7 is a lower bound of the saddle-node bifurcation curve is to consider a 
special rational parametrization of the straight line 7b — 5m = 0, m = m(s), b = b(s); 
see dSJ). With this parametrization, the coordinates of both critical points on this line 
and the eigenvalues of the saddle point are rational functions of s. These facts diminish 
the computational difficulty helping to prove the result. 

When m is large, we adapt the above approach constructing piecewise algebraic closed 
curves, which also approximate the separatrices of the saddle point and takes into 
account the region where both separatrices touch, for the first time, the negative x-axis. 
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The results that involve complicated and long algebraic manipulations are done both 
with Mathematica and Maple. 

The methods introduced in this work can also be useful to quantitatively study the 
unfolding of other singularities. For instance, the cusp of codimension 3 (our case has 
codimension 2) considered in [5], the cases considered in [HE], or the one studied by 
Takens in [TS] and quoted in [13 p. 482] could be approached with our tools. 

The paper is organized as follows. In Section [2J changes of variables in phase and 
parameter spaces are presented to shorten the computations. The new results on the 
bifurcation curve, close to the origin in the parameter space, are proved in Section [3] 
and in Section H] we prove the results for m big enough. Section [5] is devoted to prove 
Theorems [T] and [2j Finally, in Section [61 Perko's formulation of the conjectures studied 
in this paper is recalled and it is proved that they are equivalent to the ones stated 
above. 



2. Changes of variables 

As usual, to simplify the computations, we introduce changes of variables in both 
coordinates and parameters. First, we move the saddle point of system ([1]) to the origin. 
Second, we change the parameters in such a way that the eigenvalues of the saddle point 
be rational functions of the new parameters. We obtain in this way simpler expressions 
for the bifurcation curves. 

System (Tj[)) is transformed by the change of variables (x, ?/)—)■ (x + m, y) into 

X = y ' (3) 
y' = 2mx + (b + m)y + x 2 + xy. 

In order to obtain rational expressions for the eigenvalues associated to the saddle point 
we introduce new parameters M and B defined by 

M 2 = Q + m) 2 + 8m B = b + m ^ 
Then, system ([3]) writes as 

y' = (M 2 -B 2 )x + 2By + x 2 + xy. ^ 

The origin of © is a saddle point and the focus of ([I]), (— m,0), changes to (B 2 — 
M 2 ,0). 

In the parameter space the Hopf bifurcation curve b = m becomes B 2 + 2B — M 2 = 0. 
The condition of existence of an invariant straight line b = m — 1 is moved to (B + 1) 2 — 
M 2 = and the lower bound given in [9], b = m/2, writes as 3B 2 + 8B — 3M 2 = 0. 
Moreover, the origin goes to the origin and large values of m also correspond with large 
values of M. 

The homoclinic bifurcation curve b = b*(m) has a new expression for system ([3]). We 
denoteS it by W(M, B) = 0. As it is located between the curves listed in the previous 
paragraph, we have that W(M, B) = is contained in the set 

TZ = { (M, B) e M 2 : B > 0,M>0,3B 2 + 8B>3M 2 ,B 2 + 2B < M 2 , B > M - 1} , (6) 

see Figure [H 



^Although numerically it seems that we can write it as B = B*(M) we have no enough information 
on b = b*(m) to ensure this fact. 
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Without loss of generality, we consider system (j^J) only in the region 1Z. The known 
results about b = b*{m) can be translated to analogous results for system (0). 




M 



FIGURE 4. Curves that define the region 1Z (colored), that contains 
W(M,B) = 

3. Bounds near the origin 

This section is devoted to the study of lower and upper bounds for b = b* (m) for small 
values of m. The first result proves the lower bound given in Perko's Conjecture XX, 
by using a suitable rational parametrization of the straight line 7b — 5m = 0. The 
second result provides a new algebraic upper bound. In the M, B parameters, this 
upper bound, D(M,B) = is given, in an implicit way, by a polynomial D of degree 
14. When we transform it into the m, b variables, we get an algebraic curve of degree 
25. In Theorem [2] we give a much simpler upper bound, see the details in Section |5j 

Notice that the lines 7b — 5m = and b = m — 1 intersect at (m, b) = (7/2, 5/2) and 
for m > 7 we have m — 1 > 5m/7. Recall that it is known that m — 1 < b*(m), see [13J 
or Lemma O Hence, in next result, we only consider m < 7/2. 

Proposition 3. For all m < 7/2 it holds that b*(m) > 5m/7. 

Proof. Recall that system (j^J) has a saddle point at the origin and the linear approxi- 
mation of its separatrices is given by the equation 

C 2 (x, y) := (y-(B + M)x)(y - (B - M)x) = 0, (7) 

which defines two straight lines. The one with slope B + M (resp. B — M) is tangent 
at the origin to the unstable (resp. stable) separatrix. 

The curve 7b — 5m = 0, for m, b > 0, can be rationally parametrized as 

(m(s), b(s)) = ( , ^ . , , 5 ^ ^ , (8) 
V y h y 11 ^2(65 + 7) ' 2(6s + 7) J ' y J 

with s > 0. From the above expression we have that M(s) = (7s + 3s 2 ) /(6s + 7) and 
B(s) = 3s 2 /(6s + 7). Hence, with this parametrization, the slopes of the separatrices 
are rational functions of s. 

For iV G {1, 2}, let us consider an algebraic curve of degree 2(iV + 1) of the form 

2(JV+1) 

C(x,y) :=C 2 (x,y)+ ^ (c M x k + c k - hl x k ^y + c k . 2 ,2 x k ~ V) = 0, (9) 

fc=3 
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to be determined. Notice that the above curve is quadratic in y. Following the method 
described in [7j, we impose that this expression defines a curve as close as possible to 
the separatrices. It means that it should coincide at the origin with the separatrices up 
to the highest possible derivative orders. For this purpose, first we evaluate the Taylor 
series expansions of both separatrices close to the origin. We express the separatrices 
as functions of x, 

oo 

y = ^(ar) := (10) 

k=l 

where the superscript sign determines the separatrix that we approach in each case. 
Then at are real numbers obtained from the identity 



d(y-^(x)) / 7s 2 6s 2 2 

V "I o 7. « x + 7, « y + x + xy 



dx ' dy \6s + 7 6s + 7 

Straightforward computations show that the first values of are 

-7s 



= 0. 

y=*±(x) 



1 6s + 7' 

(s + l)(6s + 7) _ s-7 



2 3s (4 s + 7) ' 2 3s (2 s + 7)' 

+ _ (s + 1) (6s + 7) 2 (5s + 14) _ _ (7- s) (6s + 7) 2 (s + 2) 

° 3 ~ 18s 3 (4 s + 7) 2 (9 s + 14) ' ° 3 ~ 18s 3 (2 s + 7) 2 (3 s + 14) ' 

In fact all the coefficients of $ ± have rational expressions depending on s with non 
vanishing denominators when s > 0. 
Substituting (fit)]) in IJBjl we get 



G ± (x):=C(xM(x)) = Y J 9tx k . (11) 

fc=3 

From definition (JTj), g^ vanish for k = 1, 2. Imposing that the curve defined by C(x, y) = 
in becomes closer to the separatrices provides extra conditions = for higher 
values of k. We have 6N free coefficients Qj. We can fix these coefficients imposing 
that = 0, for k = 3,4,..., 3N + 2 and solving the system. We obtain that all the 
Ci^j are rational functions of s, which are well defined for s > 0. We call U(x, y, s) the 
numerator of C(x,y). We have 

U(x,y,s) = T 2 (x, s)y 2 + T x (x, s)y + T (x, s), (12) 

where T ,Ti and T 2 are polynomials. 

The proof continues showing that, for a given set of values of s, the algebraic curve 
U(x, y, s) — defines a positive invariant closed region that contains the focus point; 
see the right picture in Figure [31 This assertion follows if we prove: 

(I) The curve U(x, y,s) = has a loop, as it is shown in Figure El included in the 

strip x(s) < x < for a given negative value x(s). 
(II) This curve is without contact for the vector field on this strip. 
(Ill) The vector field points in on the loop. 

Notice that if we prove these properties for s G (0, 7] we ensure that the straight line 
7b — 5m = is a lower bound of b = b*(m) for m G (0, 7/2]. 

Taking the curve U(x, y,s) — corresponding to iV = 1, we can only prove the above 
assertion when s G (0,5]. When s G [5,7) we need to consider N = 2. We will detail 
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only the proof for N = 1. For the case N = 2, we will describe only the differences 
between the two cases. 

So, let us prove (I)-(III) taking N = 1 and s G (0, 5]. 

Proof of (I). Notice that the curve U(x,y, s) = can be written as 

-Ti(x,s) ± y/A(x,s) 
V ~ 2T 2 ( Xj s) 

where A := T 2 — AT 2 Tq. Straightforward computations show that 

A(x, s) = x 2 i? 4 (x, s), (13) 

where i? 4 is a polynomial of degree 4 in x and of degree 30 in s. Hence, item (I) will 
follow if we prove that there is a negative value x(s) such that: 

(i) -R 4 (x(s), s) = and R^x, s) > in (x(s), 0). 

(ii) T 2 (x,s) ^0 in [x(s),0). 

We start proving (i). For each s, the coefficients tq(s) and r 4 (s) of minimum and 
maximum degree of _R 4 in x are both positive, it has exactly two simple negative zeros 
and it has no double zeros. The non-existence of double zeros is due to the fact that 

i? 4 (s) = Res(-R 4 (x, s), <9i? 4 (x, s)/dx, x) 7^ 0, s G (0, 5], 

where Res(-,-,x) denotes the resultant with respect to x; see for instance [llj. These 
properties follow studying the roots of r , r 4 and i? 4 , that are polynomials with rational 
coefficients with respective degrees 30, 26 and 190. Their Sturm sequences ensure that 
they have no positive roots for all s G (0, 5]. We take x(s) to be the maximum of the 
negative zeros of i? 4 (x, s). 

Let us prove (ii). The polynomial T 2 (x, s) has degrees 2 and 14 in x and s, respectively. 
Straightforward computations, using its Sturm sequence, show that the discriminant 
with respect to x of T 2 is a polynomial with rational coefficients of degree 26 in s without 
positive zeros. Hence, it does not vanish when s G (0, 5]. It is easy to see that T 2 (0, s) ^ 
when s > because it is a polynomial of degree 26 with positive coefficients. For 
proving that T 2 (x(s), s) ^ for s G (0, 5] first we show that Res(i? 4 (x, s),T 2 (x, s),x) ^ 
for every s in this interval. In fact, it is a polynomial of degree 106 in s with no real 
roots in s G (0, 5]. Therefore, for these values of s, the number of real roots of T 2 (x, s) 
in [x(s),0) does not depend of s. Studying for instance the case s = 1 we can easily 
verify that T 2 (x, 1) has no real roots in \x(l), 0) and then we have the desired result. 

Proof of (II). We have to show that the vector field ([5]) is never tangent to U (x, y, s) = 
0. To prove this fact we study the common zeros between U and its derivative with 
respect to the vector field 

• dU{x,y,s) dU{x,y,s) ( 7s 2 6s 2 2 \ 

U(x,y,s):= Qx y+ ^ ^— x+— y + x + xy j . 

We get 

Res(U (x, y } s) } U (x, y, s),y) = x 12 S 4 (x, s), (14) 

where S^x, s) is a polynomial of degree 4 in a; and of degree 65 in s. We want to prove 
that S4 does not change sign for x G [x(s), 0) and s G (0, 5]. 

For a given value of s, say s — 1, the result follows directly from the Sturm method. 
After that, we proceed like in the proof of (I). We compute 



Res(i? 4 (x, s), S^x, s), x) and ~Res(Si(x,s),dS±{x,s)/dx,x). 
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We obtain two polynomials in s of degrees 362 and 438, respectively. In the interval 
(0, 5], the former polynomial does not vanish and the latter has six different real roots. 
Then, as can be seen with the Sturm method, in each of the seven intervals defined by 
these roots the relative position of the maximum negative zero of R^x, s) with respect 
to the negatives zeros of S^(x, s) does not change. Choosing s in each interval we can 
check that there are no zeros of S4 in [x(s),0) when s G (0,5], as we wanted to prove. 
Moreover, S^O, s) does not vanish on s G (0, 5]. Then the curve U(x, y, s) = is without 
contact in the region where the loop is defined. 

We remark that this step is the one that does not work in the whole interval (0, 7), 
taking N — 1, because 64(0, s) has a zero close to 5.08. 

Proof of (HI). We finish the case N — 1 checking that the vector field points in at the 
intersection point (xo(s), 0), of the loop contained in the curve U(x, y, s) — with the x- 
axis. We prove that U(x (s), 0, s), |^(x (s), 0, s) and ^(x (s), 0, s) are all negative for 
s G (0, 5], where x (s) is the biggest negative zero of U(x, 0, s). As in the above items, 
it is enough to check the inequalities for a concrete value of s. This holds because, 
straightforward computations show that Res(U(x, 0, s)/x 2 , U(x, 0, s)/x 2 , x), which is a 
polynomial of degree 69 in s, does not vanish in (0, 5], as can be verified with the Sturm 
method. 

As we have already said, the proof finishes taking N = 2 in ([9]) and following the above 
procedure for s G [5,7). All the qualitative properties remain unchanged. In fact, ex- 
pressions ffTBl and (fT4"j) change to A(x, s) = x 2 R$(x, s) and Res(U(x, y, s), U(x, y, s),y) = 
x 18 S$(x, s), respectively. Here R$ and 5*8 are now polynomials of degree 8 in x and 98 
and 229 in s, respectively. □ 

Proposition 4. For all < M < 30, the graph of W(M, B) = is below the graph of 
the function defined by the branch of the algebraic curve D(M,B) = 0, that writes as 

-383292M 14 -1910439M 13 J B-3223665M 12 J B 2 -314748M 11 J B 3 + 5603940M 10 J B 4 
+ 5541141M 9 J B 5 -2323401M 8 J B 6 -7154664M 7 J B 7 -3397092M 6 J B 8 + 2020587M 5 J B 9 
+ 3218997M 4 J B 10 + 1742052M 3 J B 11 +499644M 2 J B 12 + 76071M J B 13 +4869 J B 14 

- 500742M 13 - 787023M 12 J B +4493070M 11 J B 2 + 14795091M 10 J B 3 + 1 1566572M 9 5 4 
-11585754M 8 J B 5 -24443044M 7 5 6 -8307134M 6 5 7 + 12772706M 5 J B 8 

+ 16289545M 4 J B 9 + 8662166M 3 J B 10 + 2509195M 2 J B 11 + 388536M J B 12 + 25344 J B 13 
-174798M 12 + 1420524M 11 J B + 7005177M 10 J B 2 +4483350M 9 J B 3 -16943919M 8 J B 4 
-28501282M 7 J B 5 -3691132M 6 J B 6 + 27741570M 5 J B 7 + 31479694M 4 J B 8 + 16742014M 3 J B 9 
+4938651M 2 J B 10 + 781536M J B 11 + 52119 J B 12 -48600M 11 + 855846M 10 J B+404136M 9 J B 2 
-8473533M 8 J B 3 -14178838M 7 J B 4 + 2903273M 6 J B 5 + 26313718M 5 J B 6 + 28894211M 4 J B 7 
+ 15714446M 3 J B 8 +4769205M 2 J B 9 + 775554M5 10 + 530465 11 + 226800M 9 J B 
-1138914M 8 J B 2 -2990748M 7 J B 3 + 2116351M 6 J B 4 + 10975549M 5 J B 5 
+ 12542602M 4 J B 6 + 7156446 M 3 B 7 + 2261723M 2 5 8 + 380289MB 9 + 26766B 10 

- 264600M 7 J B 2 + 218442M 6 J B 3 + 1575182M 5 J B 4 + 2042992M 4 5 5 + 1262428M 3 J B 6 
+421554M 2 J B 7 + 73806M J B 8 + 5364 J B 9 = 0, 
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and whose series expansion at M = is 

r, 3 ~ 180,. 2366307 

B = -M M + M + 0(M 7 ). 

7 2401 90589730 v ; 

Proof. It follows using the same ideas and techniques developed in the proof of Proposi- 
tion [3] but with bigger computational difficulties. We only comment the main differences 
avoiding the details. 

Here it is sufficient to consider a curve C(x,y) = like in (EJ) of degree four. Nev- 
ertheless, its coefficients are obtained in a different way. We solve the system given by 
= 0, k — 3,4,5,6 and g% = g7 = 0. It means that we match the separatrices with 
different orders. Finally, using the same notation that in the proof of Proposition [31 the 
expression of D(M, B) = is obtained imposing the condition gj(M, B) = 0, following 
the same strategy as in [7J. □ 

Remark 5. In the m,b variables the expression of D(M, B) = provided in the latter 
result is transformed into a new algebraic curve of degree 25 with 257 monomials. 



4. Bounds up to infinity 

For the sake of completeness we also include a proof of the inequality b*(m) > m — 1, 
different to the one given in [13]. It is based on the Bendisxon-Dulac criterion, and it 
is quite simple. 

Lemma 6. For all m > 0, it holds that b*(m) > m — 1. 

Proof. It suffices to prove that, when b = m — 1, system (PQ) has no limit cycles. It writes 

as 

x' = y := P(x,y), 

y' = — m 2 + (m — l)y + x 2 + xy := Q(x, y), 

and it has the invariant straight line i := {x + y — m = 0}. Therefore the limit cycles, 
if they exist, are contained in U := M 2 \ £. We have that on U, 

\L{x,yY L(x,y) J dx\L(x,y)J dy \L(x,y)J L(x,y) 

where L(x,y) = x + y — m. Therefore, we can apply the well-known Bendixon-Dulac 
criterion ([15]) to each one of the half planes of U, obtaining that the system has no 
limit cycles, as we wanted to prove. □ 

As a first step to obtain the upper bound of b* (m) for m > 7 given in Theorem [2] we 
study the curve W(M, B) = on regions {(M,B) : M > M a } that arrive to infinity. 
The procedure is similar to the one described in Section [31 but in this case only upper 
bounds will be provided because we only have been able to obtain negatively invariant 
regions. As we will see, these bounds will be enough to prove Perko's Conjecture I. 

Suitable negatively invariant regions are much more difficult to be found that in the 
previous section. We will construct them using piecewise algebraic curves; see Figure [5j 

Proposition 7. For every real a > 0, there exists M a > ^fa > such that 
{W(M, B) = 0} C {M - 1< B < M - 1 + ^} n {M > M a }. 
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FIGURE 5. Piecewise loop for obtaining a negatively invariant region 

Proof. We already know that when B — M — 1 the system has no limit cycles, so we 
only need to prove that, for every a, there exist a curve B = M — 1 + a/M 2 and a 
value M a such that on this curve the system has a limit cycle when M > M a . So we 
will assume that B = M - 1 + a/M 2 . 

We propose a curve C(x, y) = formed by three pieces Fj of different algebraic curves 
of degree i, for i — 1,2, 3. Figure El shows the shape of the curve C(x, y) = and the 
corresponding pieces. For each i, we use the same F^ to denote the algebraic curve 
Fi(x,y) = that contains the corresponding piece. 




F\ F2 -F3 

FIGURE 6. The different pieces F^i = 1,2,3 

The first one is 

Fx := {(x, (B - M)x) : B 2 - M 2 < x < 0}. 

Notice that it is given by the segment of the straight line tangent to the stable separatrix 
that starts at the origin and ends at the point (xi, y{), where X\ = B 2 — M 2 is the x- 
coordinate of the focus. Hence Fi(x,y) = y — [B — M)x. 

We take F 2 as a portion of the quadratic curve F 2 (x, y) = 1 + ai i0 x + a Q ^y + a 2) oX 2 + 
a i,i%y + Oo,2Z/ 2 = that passes trough (xi,yi), is tangent to Fi at this point, passes 
also by the point (a; 2 ,0), with x 2 = — 3M, and coincides at this point, until second 
order derivatives, with the solution of the differential equation. More concretely, F 2 is 
the piece between (x±,yi) and (x 2 ,0). The choice of this value for x 2 is motivated in 
Remark [HJ 

The third piece F3 is contained in a cubic curve similar to the one defined in ([9]). We 
consider 

F 3 (x, y) = C 2 (x, y) + c 3fi x 3 + c 2A x 2 y + c lt2 xy 2 + c , 3 y 3 = 0. 
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This curve is tangent to both separatrices at the origin. The four coefficients of the 
homogeneous part of degree three will be fixed to get that F 3 approaches two times 
more the unstable separatrix at the origin, passes trough (x 2 ,0) and be tangent to the 
solution of the differential equation passing by this point. In fact, F 3 is the piece of the 
curve between the origin and (x 2 , 0) which is contained in the third quadrant. 

Let us study now the behavior of the vector field on the three pieces. The gradient of 

F x at (xx, yx) is (M-B, 1), the gradients of F 2 and F 3 at (x 2 , 0) are ( 3 {b+m )"^(m-b)4 ' °) 

and (-3(5 + M)(M - B)M, 0) , respectively. If M > y/a we can conclude that the 
gradient of C points to the exterior of the closed curve C(x,y) = 0. Thus, if we prove 
that the algebraic curves defined by Fi(x,y) = and Fi(x,y) = 0, the derivative of Fi 
with respect to the vector field, have no common points, we can easily conclude that 
the vector field points to the exterior of the curve C (x, y) = for M > \fa. 

The result for the segment of straight line Fx is straightforward. 

To study the vector field on F 2 we compute R 2 ■= Res(F 2 , F 2 , y). We get that 
R 2 = (x — x 2 ) 2 p4;(x, M, a), where p^ is a polynomial of degree four in x. At M = 00, this 
polynomial has an asymptotic expansion with dominant term — M~ 6 (2M + x) 3 (4M + 
x)/46656. Hence, the dominant terms of the asymptotic expansions of the corresponding 
roots are —AM and — 2M. This latter value corresponds to a triple root. A more detailed 
computation, considering the next significant term in each coefficient of p^, shows that 
the triple root splits into a couple of complex conjugated roots and a real one, Xx, with 
asymptotic expansion 

xx(M, a) = -2M + ^pM 2 / 3 + o(M 2 / 3 ). 

v 9 

Since the other real root of p^ has asymptotic expansion x 2 (M, a) = —AM + o(M), 
Xx = B 2 -M*= (M-l + ^) 2 -M 2 = -2M+l + 0(i), 

and x 2 = — 3M, it holds that for M > M 2 (a) big enough, 

x 2 (M, a) < x 2 < xx < Xx(M, a). 

Therefore, for M > M 2 (a), R 2 does not change sign on (x 2 , Xx), as we wanted to prove. 
Moreover, it is easy to see that the vector field points to the exterior of the loop on F 2 . 

The resultant R3 of F3 and F3 with respect to x is of the form y 9 qs(y, M, a), where is 
a polynomial of degree three in y. The asymptotic expansion of the ordered coefficients 
of q 3 at M = 00 shows the sign configuration [+, — , +, — ] for M > M 3 (a) large enough. 
Thus, it is clear that q 3 has no negative zeros. Then, for M > M 3 (a), R 3 does not 
vanish in the half-plane y < where F 3 is defined. Similarly to the previous cases, we 
can check that the vector field points to the exterior of C(x,y) = along F 3 . 

The proof finishes taking M a as the maximum of M 2 (a), M 3 (a) and \fa. Then we 
have the situation given in Figure [71 □ 

Next remark clarifies some details about the previous proof. Lemma [10] will follow 
directly applying the above result for a concrete value of a. The corresponding M a is 
also given. 

Remark 8. The choice of the point (x 2 , 0) = (— 3M, 0) in the proof of Proposition^ is 
motivated by some numerical computations. We wanted to choose a point on the x-axis 
that was between the first crossing point of the separatrices of the saddle point with the 
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FIGURE 7. The negatively invariant region corresponding to the piecewise loop 

negative x-axis. Let us call (P s , 0) the first crossing point for the stable separatrix, and 
(P u ,0) for the unstable one. For several values of a, using a Runge-Kutta-Fehlberg 4-5 
method, together with degree four interpolation, we obtain numerical approximations of 
the points P s and P u for every M . Figure\B(b) shows the plots of —P s /M and —P u /M 
together with the corresponding value of the abscissa of the focus point. We remark that 
for the values of a that we have checked the limit behavior is always the same. Hence, 
the asymptotic expansions at M = oo of (P s ,0) and (P u ,0) seem to be (— 2M, 0) and 
(— 4M, 0), respectively; see Figured So, our choice of x 2 = — 3M is quite natural. 





(b) 



FIGURE 8. Separatrices of the saddle point 



5. Proof of Theorems [Q and [2] 

To prove Theorem Q] we need to translate the results of Proposition [7] to the m, b 
parameters. 

Lemma 9. The change of variables (j4j) converts the curve B — M — l + a/M 2 into the 
branch of the curve 

-4m 5 - I2m% - 8m 3 b 2 + 8m 2 b 3 + 12m6 4 + 46 5 - 60m 4 - 48m 3 6 + 88mV 

+80m6 3 + 46 4 + (-192 - 16a)m 3 + (384 - 48a)m 2 b + (64 - A8a)mb 2 

-16a6 3 + (256 - 160a)m 2 - I92amb - 32ab 2 - 256am + 64a 2 = 0, 

that satisfies 



(15) 



2a a 



b = m-l + - + — 



m m z 



1 



m- 



(16) 
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when m goes to oo. Moreover, it remains below the curve b = m — 1 + 2a/m for every 
positive m and a. 

Proof. Straightforward computations give (jT5]) . To prove the last assertion it suffices to 
see that both curves do not have common points. Clearly, from (fl6|) . the result holds 
for m big enough. 

The common points of both curves are characterized by the roots of the polynomial 
obtained substituting 5 = m — 1 + 2a/m in f)15p . Removing the denominator we obtain 

q(m) := 8m 7 + (16a + 12) m 6 + (24a + 6) m 5 + (48a 2 - 12a + l) m 4 
- 8a (3a + 1) m 3 + 24a 2 (2a + 1) m 2 - 32ma 3 + 16a 4 . 

The computation of its Sturm sequence evaluated at and at +oo gives the configura- 
tions of signs [+, — , — , — , +, +, — , — ] and [+, +, +, — , — , +, — , — ], respectively, for each 
positive a. Hence, q{m) does not have positive roots, and therefore both curves have 
no common points when m > and a > 0. □ 

Proof of TheoremUi From Proposition [7] and Lemma [9j for each a > 0, the curve 
b = m — 1 + 2a /m is an upper bound of b = 6*(m), for m > m a , where m a can 
be obtained applying the transformation (j4J) to the region M > M a . Therefore, for 

m > m n , 



which implies 



m — 1 < b (m) < m — 1 H , 

m 



b*(m) — m + 1 
lim v J . — = 

m->oo 1 jm 



because a is an arbitrary positive number. □ 

Before proving Theorem [2] we need a preliminary result. 
Lemma 10. It holds that 

51 

b*(m) < m — 1 + 



20m' 

for m > m, where m pa 6.93 is the unique positive root of the polynomial 
50331648000000000 m 17 - 243269632000000000 m 16 
-2129238425600000000 m 15 - 7211878973440000000 m 14 
+111668173209600000000 m 13 + 264470739812352000000 m 12 
-130466347912396800000 m 11 - 9197101546824499200000 m 10 
-63029001 12535388160000 m 9 + 6325778059290335232000 m 8 
+2289016716587559936000 m 7 - 46572462911915224012800 m 6 
+8515659923453703340800 m 5 - 4901243812728523876800 m 4 
-45716337137659722706080 m 3 + 6052551315638078774880 m 2 
-8203038242422388605200 m - 7953608649353382254007. 

Proof. Proposition [3 with a = 51/40, ensures that 
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defines an upper bound for W(M, B) = when M > M51/4Q. Particularizing the analysis 
done in the proof of this proposition to this particular a we can obtain an explicit value 
^5i/40- It can De taken as the largest positive root of the equation 



196608000000M 17 - 1441792000000M 16 - 535756800000M 15 
+1480294400000M 14 - 1310515200000M 13 - 1151979520000M 12 
+ 1314478080000M 11 - 727741440000M 10 - 273666816000M 9 
+443460096000M 8 - 458441856000M 7 + 61550064000M 6 + 227310753600M 5 
-162364824000M 4 - 41403030120M 3 + 82806060240M 2 - 17596287801 = 0, 



(17) 



which is approximately 7.58. 

Then the result follows from Lemma [9J Moreover, the polynomial that defines m is 
obtained computing the resultant with respect to b of (IT5"|) and the polynomial corre- 
sponding to (fT7|) translated to the variables m and b. □ 

It is useful to introduce the following notation for the lower and upper bounds of 
b*(n) given in Theorem |2] 

, / , /5m A , i / \ ■ f(5 + %m)m 25 \ /10 . 

Mm) := max m — 1 and b u (m) := mm ^4= , m — 1 + - — . (18) 

V ' V 7 / V 7+ gm ' 7my v ' 

Notice that the non-differentiability points of be and b u are at m = 7/2 and m = 7, 
respectively. 

Proof of Theorem^ Proposition [3] provides the lower bound given in the statement, 
when m < 7/2. For m > 7/2, Perko gives a proof in [T3]. A different one is presented 
in Lemma |HJ 

The proof of the second part is done by comparison of the curves in the statement 
with the ones provided by Proposition H] and Lemma [TUl 

For m > 7, b = b u {m) is an upper bound from Lemma [10] because in < 7 and 
25/7 > 51/20. 

When m < 7 the proof starts translating the curve D(M,B) = given in Proposi- 
tion H] to a new algebraic curve E(m, b) = 0, of degree 25 with 257 monomials. Now we 
compare the curves (84 + 37m)6 — (60 + 37m)m = 0, corresponding to b u {m) = 0, and 
E{m, b) = when m, b > 0. The resultant with respect to b of both polynomials takes 
the form m 15 p 34 (m), where P34 is a polynomial of degree 34 with a unique positive zero, 
mi, as can be easily seen from the Sturm method. Hence, the curves only intersect at 
(0, 0) and (mi, 61), where mi ~ 7.1 A local study of the curves close to the origin shows 
that b = b u (m) is above E(m,b) = 0. Hence, as the relative position of the graphs of 
both curves does not change when < m < 7, b = b u {m) is also an upper bound of 
b = b*(m) in the full interval. □ 

Corollary 11. Setb(m) = (b u (m) +bt(m))/2, where be andb u are given in (I18p . Then, 
the absolute and relative errors when we approximate b*{m) by b(m), are 



~ 37 

max 16* (m) — b(m)\ < < 0.31 and max 

m>o 1 v 1 v n 122 m>0 



b*(m) — b(m) 



b*(m) 



< — <0.13. 
305 
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Proof. It is not difficult to see that the maxima of the functions b u — be and (b u — be) /be 
are both at m = 7/2. Then 



max I b* (m) 

m>0 



b(m)\ < max 

m>0 



m 



6 u (7/2) - 6,(7/2) 



and 



max 

m>0 



b(m) 



< max 

m>0 



m 



bAm) 



2b e (m) 



b u {7/2)-b e {7/2) 



2be(7/2) 



37 
122 



37 
305' 



□ 



b*(m) 

Hence, the corollary follows. 

Remark 12. From Theorem^ the family of rational functions 

, 5 + (3m 
7 + [3m 

for every (3 > 0, approximates b*(m) in neighborhoods of the origin and the infinity 
simultaneously. Nevertheless, the upper bound given in Theorem [H, that corresponds 
to (3 = 37/12, is changed in a neighborhood of infinity because the family of functions 
b = m— 1 + 7/777, is a much better approximation for m big enough. The concrete values 
of (3 and 7 are fixed imposing the continuity ofb u and searching nice expressions for the 
statement of Theorem^ These values could be changed to obtain slightly better upper 
bounds. 

6. Original formulation of Perko's Conjectures 
In [131 E], Perko considers the following expression of the Bogdanov-Takens system 

(19) 

1) + fiiV + fi 2 UV, 

and proves that it has a homoclinic saddle-loop if and only if \i\ = h(fi2) for some odd 
analytic function h. Since h(—/j,2) = —h(fi 2 ) it suffices to study h either for /i 2 > or 
for jj,2 < 0. 

Original formulation of Perko's Conjectures (|13j): Let /xi = h(fx 2 ) the function 
that gives the saddle- loop bifurcation curve for system (JT9"j) . Then, 

(X) the curve /ii = h({i 2 ) is asymptotic to the hyperbola \i\ = —1/^2 for large |/i 2 |, 
i.e., fi 2 h(fi2) + 1 = 0(l//i 2 ), as /i 2 -> 00, 
(XX) for each ji 2 < 0, < /i(/i 2 ) < min{— /i 2 /7, — l//i 2 }. 

To see that the above conjectures are equivalent to the ones stated in the introduction 
we will transform system ([T]) into f[T9~|) . If we apply the change of variables u = (x + 
m)/(2m), v = y/(2m) 3 / 2 and consider the new time s = \/2mt, system ([1]) writes as 



u 
v' 



u(u 



u 



u(u - 1) + 



m 



/ 2m 



-v + \/2muv. 



Hence, we have the following equivalence among the parameters m,b and /ii,/i 2 : 

b — m , 

111 = , (jl 2 = V2m, 



\/2m 

and every curve of the form /ii = /(/i 2 ) is transformed in the variables m and b into 
b = m + y/2m /(\/2m). In particular, the curves /ii = 0, /ii/x 2 = —1 and /ii = — /x 2 /7 
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are transformed into the straight lines b = m,b = m — 1 and b = 5m/7, respectively. 
This fact shows that both Conjectures XX are equivalent. 

To compare both Conjectures X, notice that //1//2 = b — m. Therefore, we can write 
/X2A*i + 1 = 0(1/112) as b — m + 1 = 0(l/y/m), as we wanted to prove. In fact, notice 
that 

b*{m) =m + y/2171 h{\/lm). 

Indeed, since h is analytic and odd this equality proves that b*(m) is analytic in m. 

It is worth to mention that there is a third conjecture in Perko's work: For ^2 < 
the function h(n2) has a unique maximum. The tools introduced in this paper seem not 
to be adequate to approach this question. 
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